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ABSTRACT 

Starting from a complete sample of type I AGN observed by INTEGRAL in the 20-40 
keV band, we have selected a set of 8 AGN which can be classified as radio loud objects 
according to their 1 .4 GHz power density, radio to hard X-ray flux flux density ratio and radio 
morphology. The sample contains 6 Broad Line Radio Galaxies and 2 candidate ones. Most of 
the objects in our sample display a double lobe morphology, both on small and large scales. 
For all the objects, we present broad-band (1-110 keV) spectral analysis using INTEGRAL 
observations together with archival XMM-Newton, Chandra, Swift/XRT and Swift/BAT data. 
We constrain the primary continuum (photon index and cut-ofF energy), intrinsic absorption 
and reprocessing features (iron line and reflection) in most of the objects. The sources anal- 
ysed here show remarkable similarities to radio quiet type I AGN with respect to most of the 
parameters analysed; we only find marginal evidence for weaker reprocessing features in our 
objects compared to their radio quiet counterparts. Similarly we do not find any correlation 
between the spectral parameters studied and the source core dominance or radio to 20-100 
keV flux density ratios, suggesting that what makes our objects radio loud has no effect on 
their high energy characteristics. 
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1 INTRODUCTION 

Most of what is known about Active Galactic Nuclei (AGN) is es- 
sentially based on studies of radio quiet sources, which make up 
almost 90% of the entire AGN population. It is now widely ac- 
cepted that active galaxies are powered by accretion onto a su- 
permassive black hole: the observed radiation, spanning the en- 
tire electromagnetic spectrum, is produced by a cold accretion disk 
and by a hot corona, as pro posed in the so-called two-phase model 
jMaraschi & Haardt|[T997b . In the X-ray domain, the emission of 
radio quiet Broad Line AGN (Seyfert 1) is, to the first order, well 
described by a power law of photon index 1.8-2.0, extending from 
a few keV to over 100 keV; at higher energies there is evidence for 
an exponential cut - off, the exact value of which is still uncertain 
jPerola et alj|2002l : lRisalitill2002h . Secondary features such as the 
Fe Ka line and the Compton reflection component are also com- 
monly observed; they are considered to be the effects of reprocess- 
ing of the primary con tinuum and are relatively well understood 
jMushotzkv et aljll993l) . Studies performe d on samples of Broad 
Line Radio Galaxi es (BLRG from now on; ISambruna et al. I I 19991 : 
iGrandi et alj2006h have shown that these objects show optical/UV 
continuum and emission line characteristics similar to their radio 



quiet counterparts, but display some fundamental differences in 
their X-ray behaviour. BLRG, in fact, exhibit flatter/harder power 
law slopes than radio quiet Seyfert 1 galaxi es and are also known 
to have weaker reproc essing features (e.g. ISambruna et"ai]|l999l 
an d iGrandi et all2006l) . The origin of these differences is however 
far from being understood. A possible cause for the observed di- 
versity could be ascribed to a different disk geometry and/or ac- 
cretion flow efficiencies, or to the presence of jets and beaming 
effects that contaminate/dilute the AGN component and the re- 
processing features; ionised reflection, which naturally produces 
weake r reflection features, has also b een considered in the liter- 
ature teallantyne. Ross & Fabianll2002 ). Now that a large sample 
of AGN detected above 20 keV by INTEGRAL is available, it is 
possible and important to perform a comparison between different 
classes of AGN, by measuring the shape of the primary continuum 
together with the high energy cut-off and the reflection component. 
In the present work we focus on the broad-band (1-110 keV) spec- 
tral analysis of a sample of eight radio loud type 1 AGN, com- 
bining XMM-Newton, Chandra, Swift/XRT and BAT data together 
with INTEGRAL/ISGRI measurements. We also compare our re- 
sults to a sample of radio quiet Seyfert Is detected by INTEGRAL 
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and recently studied bv lPanessa et alj flOOSh and finally discuss the 
implications of our findings. 



2 THE SAMPLE 

The definition of radio loudness is rather vague, with difi^erent 
criteria applied in the literature (see IZamfir. Sulentic & Marzianil 
[2OO8 for a critical discussion). Traditionally, empirical bound- 
aries were set by the radio power density at 5 GHz (Psghz; 
iMiller. Peacock & MeadI [T990I) and the radio ( at the same fre- 
quenc y) to optical B band flux density ratio (Rn; iKellerman et"ail 
Il989() : according to these boundaries, an AGN is radio loud if 
P5G//;>10^^^ erg s"' Hz"' and Log(RB)^l. In time, different surro- 
gate definitions of both parameters have appeared in the literature 
involving radio data at different frequencies and comparison with 
even UV and X-ray fluxes to deal w ith objects selected in bands 
other than the optical. More recently [Zamfir. Sulentic & Marzianil 
( l2008h suggested the use of the source radio morphology as a fur- 
ther criterion to divide radio loud from radio quiet objects. Since 
FR I morphology is quite rare in broad line radio loud AGN, while 
double-lobe FR II appearance is much more common, these authors 
suggest the adoption of this last morphology as a further discrim- 
inator, i.e. FR Il-like objects are by default radio loud. We there- 
fore selected our radio loud AGN on the basis of similar criteria: 
the power density at 1.4 GHz (Pi 40//;), the radio (at the same fre- 
quency) to the 20-100 keV flux density ratio (R/^x) and the source 
morphology, again at 1.4 GHz. Sources were selected from a set 
of 35 Broad Line (BL) galaxies extracted from a complete sample 
of INTEGRAL detected AGlsQ Because many of these objects are 
newly discovered, it is difficult to have information for all of them 
at 5 GHz and/or in the B band, hence the use of the more ready 
available 1.4 GHz and 20-100 keV data. 1.4 GHz information has 
been obtained from the HEASARC databas^B for 29 objects in the 
sample, except in the case of IGR J 13 109-5552, where the 1.4 GHz 
flux was extrapolated using the 4.85 and 0.83 GHz fluxes avail- 
able in the HEASARC radio catalogues. INTEGRAL data were de- 
rived from the fluxes reported in Bird et al. (2007; see also table 
lA) and are available for all objects in the sample. Candidate ra- 
dio loud objects were selected on the basis of their location in the 
diagram shown in figure [T] where the radio quiet/loud boundaries 
on both axes are set by the weakest FR II source in our sample, 
i.e. IGR J2I247-I-5 058. This source has been fully discussed by 
iMolina et alj ( |2007|) and it is now a well established case of a radio 
loud AGN on the basis of various ar guments. In particula r it has 
been classified as an FR II galaxy by iMolina et al. I ( l2007h with a 
L5c«:~10" erg s^'Hz"', an order of magnitude lower than the tra- 
ditional limit, making this source one of the weakest FR II galaxies 
reported to date. 

By means of these boundaries we identified 8 radio loud AGN: 
QSO B0241-I-62, B3 0309+411, 3C 111, 3C 390.3, IGR J13109- 
5552, 4C 74.26, S5 2116+81 and IGR J21247+5058 (see table 1); 
note the ambiguous location of two additional objects, Markarian 
6 and MCG+08-1 1-1 1, which are below the adopted boundaries 
but very close to IGR J21247+5058. Information at 5 GHz and in 



* Th i s sample extracted from the 3"* IBIS/ISGRI catalogue jBird et al] 
l2007t iMalizia et alj|2008h contains all AGN detected with a significance 
greater than 5.5o" in the 20-40 keV band and having z<0. 14; it contains 
35 type 1 objects, 30 type 2 sources, 3 Narrow Line Seyfert Is, and 7 
Blazai-s/QSOs 

^ http://heasarc.gsfc.nasa.gov/ 



the B band are available for all these sources so that a cross-check 
can be made on their radio power density and radio-to-optical flux 
density ratio acco rding to the more conventional definition used by 
iMiller. Peacock & Mead ( 1990) and Kellerman et al. ( 1989). Both 
values are listed in table 1 for all 8 candidates plus Markarian 6 and 
MCG+08-1 1-01 1, which we keep under consideration due to their 
location in figure [T] 5 AGN have Pjo/zi^lO^^ erg s"' Hz"' and in 
any case all but Markarian 6 and MCG+08-1 1-01 1 have a 5 GHz 
power density exceeding that of IGR J21247+5058. All the objects, 
including Mrk 6 and MCG+08-1 1-011, have Log (Rb)> 1 ; note that 
in IGR J21247+5058, Rb cannot be estimated due to the lack of 
information on the B magnitude. 

Next we examined the radio morphology of all our 8 radio 
loud candidate AGN, plus again Mrk 6 and MCG+08-1 1-01 1 (see 
table 1). All but 2 sources (i.e. QSO B0241+62 and IGR J13I09- 
5552) are well known Broad Line Radio Galaxies (see table 1 and 
NED, Ballantyne 2005 and Molina et al. 2007) and show a mor- 
phology with two lobes extending from the central nucleus (see 
for example NVSS contour maps available in NED); the extension 
of the lobes for B3 0309+411 and 4C 74.26 is m ore than 1 Mpc 
hence their classification as giant radio galaxies ide Bruvrilll989l ; 
[kiley et al. 1988). These double-lobed objects are classified in the 
literature as FR II galaxies, with the exception of 4C 74.26 and 
S5 2116+81, which still have an uncertain/unknown nature. QSO 
B0241+62 and IGR J13109-5552 have been poorly studied at radio 
frequencies and their classification is uncertain; both are Seyfert 1 
according to NED and Masetti et al. 2008. QSO B0241+62, shows 
a compact unresolved core in the NVSS map, although high reso- 
lution radio imaging indic ates a possible double- lobe morphology 
within this core structure j Tzanetakis et aljri978l) . The sky region 
containing IGR J13109-5552 has not been mapped by the NVSS, 
but the source is present in the Molonglo Galactic Plane Survey 
2nd Epoch (MGPS-2) Compact Source Catalogue jMurphy et al.l 
2007); a cut-off image of the source shows an unresolved but elon- 
gated structure. QSO B0241+62 is reported as a relatively strong 
source i n various radio surve ys, has a flat radio spectrum with a 
Q'~O.10 jVollmer et al.ll20()^ and so almost certainly qualifies as 
a radio loud AGN. IGR J13 109-5552 has less information at radio 
frequencies but the few available data provide a flat spectrum also 
in this case (q'~0.3). 

We also checked, for completeness, the morphology of 
Markarian 6 and MCG+08-11-011 in the NVSS/NED database: 
both are similar to QSO B0241+62, i.e. display a compact 
core which, at higher resolution, reveals a clear (in Markar- 
2006) or possible (in MCG+08-11-011; 



ian 6; Kharb et al. 



lUlvestad & Wilsonlll9 86) double l obe structure. Howe ver their ra- 
dio spectra are steeper (g~0.6-0.7; IVollmer et al.ll2005h than in the 
previous two cases and more similar to radio quiet AGN. This evi- 
dence, coupled with their lower radio power argues against a classi- 
fication as radio loud AGN and therefore both objects are not con- 
sidered further. 

Table lA lists all relevant information on the 8 sources 
chosen for this study, including redshift. Galactic column den- 
sity in their direction, 20-100 keV flux and luminosity and the 
5/8 GHz radio core d ominance (CD). The Galac tic column den- 
sities are taken from iDickey & LockmanI 
rameter CD=Score/(Sioi-Score) is taken from 



199lh. while the pa 



Fan & Zhan j ( l2003h ; 



the only exceptions are B3 0309+411 and IGR J21247+5058, 
for which the CD parameter has been evaluated using data from 



^ SyOCV 
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Figure 1. 1.4 GHz power density versus the 1.4 GHz to 20-100 keV flux density ratio for most sources in the INTEGRAL BL AGN catalogue. The boundaries 
on both axis are set by the values of the weakest FRII in the sample (IGR J21247+5058) and are defined by the lines drawn in the figure. 



Ciindidate Radio Loud AGN 



Name 


Type 


L«8fP5G//,7* 


Log(Re) 






erg s~ ' Hz~ ' 




QSO B0241+62 




31.230 


2.8-2.9 


B3 0309+411 


BLRG. FRII? 


32.322 


3.3 


3C 111 


BLRG. FRII 


31.919 


4.4-5.0 


IGR J13109-5552 




31.929 


2.7 


3C 390.3 


BLRG. FRII 


32.477 


3.7 


4C 74.26 


BLRG. FRI/FRII 


31.939 


1.9-2.0 


S5 2116+81 


BLRG. FR? 


31.771 


1.9 


IGR J2 1247+5058 


BLRG, FRII 


31.041 




MCG+08-11-01I 




29.886 


1.6 


Mrkfi 




29.934 


1.5-1.8 



llshwara-Chandra ~&. Saikij ( Il999h and lMolina et alj j2007h respec- 
tively. No information on CD is available for IGR J13109-5552. 

The radio core dominance spans a large range of values from 
0.04 to 1.9. Sources with high values of this parameter are those 
in which the beamed radio emission is emitted in a direction closer 
to the line of sight; in these sources, jet emission is likely to play 
a greater role than in sources with smaller values of CD. Contrary 
to expectations, we do not find any convincing correlation between 
CD and either (P\.aghz) or ^hx, indicating that in at least some of 
our sources a jet component alone cannot be the origin of the source 
radio loudness. As a final remark, we point out that our 8 radio loud 
sources make up around 20% of the 20-40 keV complete sample of 
BL active galaxies and are therefore a non-negligible fraction of the 
INTEGRAL AGN population. 



3 DATA REDUCTION 

The sources in our sample have all been observed by XMM-Newton 
with the exceptions of QSO B0241+62, for which Chandra data 
are available, and IGR J13109-5552 and S5 2116+ 81, for which 
Swift/XRT data have inste ad been used. MOS and pn JStriider et al.l 
I2OO1I ; iTurner et al.ll200l[) data were reprocessed using the XMM- 
Newton Standard Analysis Software (SAS) version 7.0 employing 
the latest available calibration files. Only patterns corresponding 
to single, double, triple and quadruple X-ray events for the two 
MOS cameras were selected (PATTERNS 12), while for the pn only 
single and double events (PATTERN<4) were taken into account; 
the standard selection filter FLAG=0 was applied. Observations 
have been filtered for periods of high background and the result- 
ing exposures are listed in table 2. Source counts were extracted 



from circular regions of typically 40"-50" of radius centered on 
the source, while background spectra were extracted from circular 
regions close to the source or from source-free regions of typically 
20" radius. The ancillary response matrices (ARFs) and the detec- 
tor response matrices (RMFs) were generated using the XMM-SAS 
tasks arfgen and rmfgen; spectral channels were rebinned in order 
to achieve a minimum of 20 counts per each bin. For sources af- 
fected by pile-up (3C ill, 4C 74.26 and 3C 390.3) the central 5" of 
the PSF have been excised and spectra have been extracted from an- 
nular regions of typically 50" external radius. XRT data reduction 
for IGR J13109-5552 and S5 2116+81 was performed using the 
XRTDAS vl.8.0 standard data pipeline package (xrtpipeline v. 
0.10.3) in order to produce screened event files. All data w ere col- 
lected in the Photon Counting (PC) mode jHill et al.ll2004l) . adopt- 
ing the standard grade filtering (0-12 for PC) according to the XRT 
nomenclature. Source data have been extracted using photons in 
a circular region of radius 20"; background data have been taken 
from various uncontaminated regions near the X-ray source, using 
either a circular region of difi"erent radius or an annulus surrounding 
the source. 

Chandra data reduction for QSO B0241+62 was performed 
with CIAO 3.4 and CALDB 3.2.4 to apply the latest gain correc- 
tions. Subsequent filtering on event grade and exclusion of periods 
with high background resulted in a total exposure of 34 ks. The 
CIAO script "psextract" was used to generate the source spectrum, 
with its appropriate background and response files; spectral data 
were extracted from a circular region of radius ~ 4" while back- 
ground files were generated using a region of -18" in diameter 
With a count rate of 0.14 counts/frame, the pile-up fraction is in- 
significant (<10%). 

The INTEGRAL data reported here consist of several pointings 
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Table lA 



Source Data 



Name 


RA 


Dec 


z 


10--cm-2 


Ifl-l^erg cm"-s"l 


^2a-\mkeV 
lO'^'^erg s"' 


CD- 


QSO B0241+62 


41.285 


+62.480 


0.044 


0.75 


6.34 


2.63 


1.94 


B3 0309+41 1 


48.273 


+41.343 


0.136 


0.13 


2.57 


9.75 


1.17 


3C 111 


64.573 


+38.014 


0.0485 


0.32 


9.46 


4.75 


0.04 


IGRJ13109-5552V 


197.682 


-55.863 


0.085 


0.22 


2.42 


0.33 




3C 390.3 


280.586 


+79.781 


0.0561 


0.04 


6.06 


4.06 


0.09 


4C 74.26 


310.585 


+75.145 


0.104 


0.12 


4.87 


10.97 


0.95 


S5 2116+81 


318.492 


+82.072 


0.084 


0.07 


4.06 


6.02 


1.78 


IGRJ21247+5058 


321.172 


+50.972 


0.02 


1.11 


10.90 


0.94 


0.76 



^ : radio core dominance at 5 GHz, except for B3 0309+41 IB for which it is measured at 8 GHz (see references in the text). 



perfo rmed by the IBIS/ISGRI l lUbertini et all l2003l : iLebrun et al.l 
l2003h instrument between revolutions 12 and 429, i.e. the period 
from launch to the end of April 2006. ISGRI images for each avail- 
able pointing were generated in various energy bands using the 
ISDC offline scientific analysis software OSA jGoldwurm et al.l 
l2003h version 5.1. Count rates at the position of the source were 
extracted from individual images in order to provide light curves 
in various energy bands; from these light curves, average fluxes 
were then extracted and combined to produce an average source 
spectrum (see iBird et al. 2007. for details). Analysis was performed 
in the 20-110 keV band. In the present analysis we also made 
use of public Swift/BAT spectra, retrieved on the wet|3; spectra are 
from the first 9 months of operations of the Swift/BAT telescope 
teaumgartner et al.ll2 008). 

Table 2 reports the observation log of each source, i.e. the ob- 
servation date of each X-ray observation and the exposure time for 
each instrument employed in the spectral analysis, with the excep- 
tion of Swi/r/SAJ. 



4 BROAD-BAND SPECTRAL ANALYSIS 

The XMM-Newlon, Chandra and Swift/XRT data were fitted to- 
gether w ith INTE GRAL/ISGRI and Swift/BAT data using XSPEC 
v.11.2.3 ( lAmaud[l996 ): errors are quoted at 90% confidence level 
for one parameter of interest (A;jf^=2.71). In the fitting procedure, 
a multiplicative constant, C, has been introduced to take into ac- 
count possible cross-calibration mismatches between the X-ray 
and the soft gamma-ray data. When treating the INTEGRAL data, 
this constant Ci has been found to be close to 1 with respect 
to XMM-N ewton, Swift/XRT and Chandra using various source 
typologies fcandi et al.1 12007| : iDe Rosa et"^ |2008| ; iMasetti et al] 
I2OO7I : IPanessa et alj|2008l) . so that significant deviation from this 
value can be confidently ascribed to source flux variability; also 
when considering BAT data, a cross-calibration constant C2 differ- 
ent to 1 often means flux varia tion in the source a nalysed (see for 
example IXiello etal . 2008 and lWinter et alj|2008h . We also intro- 
duced cross calibration constants between the XMM instruments, 
pn/MOSl and MOS2/MOS1; these were left free to vary and al- 
ways found to be in the range 0.97-1.04. In all our fits. Galactic 
absorption (see table lA) has already been taken into account so 
that any column densities reported in subsequent tables refer to ab- 
sorption intrinsic to the source. 



fittpV/swift. gsfc.nasa.gov/docs/swift/results/bs9mon/ 



4.1 Towards a first approximation of tlie broad-band 
continuum. 

We firstly fitted the broad-band data in the 0.5-110 keV energy 
range employing a simple power law (absorbed only by Galactic 
column density) in order to identify typical features of the AGN 
spectra. Some of our objects show evidence of excess counts be- 
low 1 keV; however, since the study of the soft excess is not the 
main objective of the present work, but rather the understanding of 
the high energy emission characteristics (photon index, high energy 
cut-off and reflection), we have focussed the analysis in the 1-110 
keV energy range. Starting from the simple power law we intro- 
duced various spectral features such as intrinsic absorption and iron 
line, each time performing an F-test to verify the statistical signifi- 
cance of each new component and to provide a basic model able to 
describe the data at least to a first approximation (see table 3 and 
residuals with respect to this model in figures (2] [3]|1] and O. 

Inspection of figure|4](left panel) suggests that the model em- 
ployed in table 3 does not provide a good description of the data 
of 3C 390.3, where the residuals have an unusual concave shape, 
which may be due to the presence of a low energy component (soft 
excess) still affecting the spectrum around 2-3 keV. To take this into 
account we have added an extra feature to the model in the form of 
another power law or a blackbody component: both fits provide 
an equally significant improvement (typically more than 99%) and 
more acceptable residuals. The blackbody temperat ure is slightly 
highe r than typically observed in Seyfert 1 galaxies dPanessa et alj 
'2008!). It is interesting to note that so far no soft excess has been re- 
ported for this source exce pt for an old claim never confirmed after- 
wards i lGrandietal.ll 19991) : however a Chandra image of 3C 390.3 
reveals extended soft X-ray emission around the nucleus so the 
soft excess photon s we detect co uld be related to this compo nent 
dKadler et alj2004 ' and figure 7 in lHardcastle & Crosto n 20^. We 
also point out that this source is known to have a variable column 
density, as found by Grandi et al. (1999); however in our analysis 
we were not able to detect any absorption in excess to the Galac- 
tic one. Due to the evidence found, in the following we adopt the 
model wag*wa*(bb+po+zga) (table 3A and figure |6]right panel) 
as our basic representation of 3C 390.3. 

It is also evident from table 3 and figure [5] (right panel) that 
IGR J21247+5058, has a very flat power law slope and conse- 
quently shows a mismatch between X-ray and soft gamma-ray data. 
It is however known from the literature that this so urce requires 
a model that takes into account complex absorption iMolina et al.l 
|2007) in the form of two layers of material partially covering the 
central source. Hence, in the case of IGR J21247+5058, the sim- 
ple intrinsic photoelectric absorption has been substituted by two 
partial covering components (wag*pcfabs*pcfabs*(po+zga) in 
XSPEC, see results in table 3B); when this model is applied, the 
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Table 2 



Observations Log 



Name 


Obs. Date 


Exposure (pn) 

(sec) 


Exposure (MOS) 
(sec) 


Exposure (XRT) 
(sec) 


Exposure {Chandra) 
(sec) 


Exposure (INTEGRAL) 
(sec) 


QSO B0241+62 


06-07/04/2001 








34000 


371986 


B3 0309+411 


04/09/2005 


10990.1 


19703.1/16138.6 






506677 


3C 111 


14/03/2001 


14542.1 


76262.3/7116.3 






200566 


IGR J13109-5552 


25/12/2006 






6.58 




1288 


3C 390.3 


08/10/2004 


38762.3 


32450.7/34603.9 






305807 


4C 74.26 


06/02/2004 


25643.5 


28960.4/28960.4 






121316 


S5 2116+81 


25/06/2006 






51114 




198462 




17/10/2006 






52374 






IGR J21247+5058 


06/1 1/2005 


22166 


24245/23614 






1067000 



: note thai for INTEGRAL the expo.sure refer to a number of pointings in the period between launch and April 2006. 



power law index becomes steeper and the X-ray/gamma-ray match 
is more acceptable (see figure|6]left panel). 

Using the results reported in table 3, 3 A and 3B we can al- 
ready draw some conclusions, as the models employed at this stage 
provide, to a first approximation, a good description of the broad- 
band spectra of all our sources. In all but 2 AGN (B3 0309+41 1 and 
3C 390.3), an intrinsic absorption component is strongly required 
by the data, at more than 90% confidence level for IGR J13109- 
5552 and S5 2116+81 and at more than 99% confidence level for 
the remaining objects. Note, however, that in 3C 111 the measured 
absorption coul d be due to a mol ecular cloud located between us 
and the source dUngeretalJI 19851) and not intrinsic to the AGN. It 
is also interesting to note that only one out of eight objects displays 
complex absorption; a similar fraction was found in a sample of 9 
radio quiet Seyfert 1 galaxies discussed by IPanessa et alj ( 1200 Sl) . 
where at least two objects required one or two layers of cold mate- 
rial partially obscuring the central nucleus. 

The Fe Ka line is detected in all the sources (99% confidence 
level), except in IGR J13109-5552 and S5 2116+81, but this may 
be due to XRT not being sensitive enough to detect the cold iron 
line around 6.4 keV. The line is narrow in all of our sources (hence 
the line width has been fixed to 10 eV), except in 4C 74.26, in 
whi ch a quite strong a nd p ossibly broad iron l ine is detected (see 
also lBallantvnell2005l and iLarsson et al.ll2008l) . The line equiva- 
lent widths (EW) are found to be <100 eV in all of our AGN. In 
two sources, namely 3C 1 1 1 and IGR J21247+5058, the equivalent 
width values are close to the capability limits of moderate reso- 
lution CCD instruments like XMM-Newton, calling for some cau- 
tion in distinguishing between a real feature or just local noise. If 
the line is related to reflection, the observed values would imply 
a low reflection component parameter {K=Q.I2t[): in this case one 
may expect to find the EW to R ratio of the order of 100-130 eV 
jPerola et alj2002h and consequently values of R^^i. We will come 
back to this point later. 

Finally a note on the cross-calibration constants between the 
X-ray data and IBIS/BAT points: in some sources these constants 
are not consistent with 1, indicating, as expected, flux variations 
between the X-ray snap-shot observations and the time-averaged 
high energy spectral data. 



around 100 keV for 3C 1 1 1 and 80 keV for IGR J21247+5058. As 
far as the other sources are concerned, we could only set a lower 
limit for their cut-off energies at around 60-80 keV, apart from 3C 
390.3, where the high energy cut-off could not be constrained at 
all. It is also interesting to note that the values of Ci and C2 tend to 
be higher when a cut-off energy component is added to the simple 
power law. 

4.3 Incidence of the reflection component. 

We then introduced reflection in our basic model of tables 3, 3A 
and 3B to verify the incidence of this spectral component in our 
broad-band data (again refer to tables 5, 5 A and 5B for model de- 
scription in each source case). We substituted the power law with 
the pexrav model in XSPEC, fixing the inclination angle at 30° (i.e. 
a nearly face-on geometry as expected in type 1 AGN) and the cut- 
off energy at 10000 keV. Taking into consideration the fact that 
the iron line EW is always <100 eV, we assume that the reflection 
parameter R cannot have very high values; for this reason we ini- 
tially allowed the reflection component to vary in the range 0-2. 
B3 0309+411 and 4C 74.26 require reflection at more than 99% 
confidence level, while 3C 390.3 and QSO B0241+62 at a lower 
confidence level (98% and 92% respectively). The reflection com- 
ponent is also significantly required (at more than 99% confidence 
level) in 3C 111 and IGR J21247+5058, the two AGN which also 
have a high energy cut-off; since in these two objects the iron line is 
very weak or possibly even absent we expect to observe negligible 
or low reflection and we thus suspect that a better fit is achieved 
when both R and E„,, are left free to vary (see next section). In S5 
2116+81 and IGR J13 109-5552 the data only provide a loose lower 
limit on R, likely due to the low statistical quality of the X-ray data. 

The reflection values are all above 0. 1 and, when constrained, 
are typically around 1 . Despite being highly required, R is uncon- 
strained in B3 0309+411 due to the upper bound of 2 set on the 
reflection; when this bound is removed, R is found to be 2.44^| 
and the constant becomes C2=2.30^q"55- It is also important to note 
that R tends to be higher than predicted from the measured iron 
line EW and that the pexrav model provides, as expected, smaller 
values of C\ and C2 with respect to the previous two models. 



4.2 Incidence of tlie high energy cut-off. 

We next introduced a high energy cut-off in the primary power law 
of our basic model (see tables 4, 4A and 4B for model description 
in each source case). Only two of the eight sources in the sample 
strongly ( at more than 99% confidence level) require a cut-off en- 
ergy, i.e. 3C 111 and IGR J21247+5058; the cut-off is localised at 



4.4 Constraining both reflection and cut-off. 

Finally, the broad-band spectra have been fitted with a cut-off 
power law reflected from neutral material (pexrav model but with 
E„„ free to vary, see table 6 and figures |7] [8] and |9] for model de- 
scription in each source case). We removed IGR J13 109-5552 and 
S5 2116+81 from the sample, as in both sources the quality of the 
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QSO B0241 +62 



B3 0309+41 1 B 




channel energy (keV) 
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Figure 2. Model to data ratios for QSO B0241+62 (left panel) and B3 0309+41 IB (right panel). The model employed is a simple power law absorbed by both 
Galactic and intrinsic column density (only in the case of QSO B0241+62) plus a narrow Gaussian line component (see Table 3). 



IGR J13109-5552 



channel energy (keV) 



channel energy (keV) 



Figure 3. Model to data ratios for 3C 111 (left panel) and IGR J13109-5552 (right panel). The model employed is a simple power law absorbed by both 
Galactic and intrinsic column density plus, in the case of 3C1 1 1, a narrow Gaussian line component (see Table 3). 



data allows no more complex fits than a simple power law (as in 
figure [Tot. Because in B3 0309+411 the reflection parameter ex- 
ceeds the value of 2, we raised the upper bound on R. Although this 
model does not always provide a significant improvement with re- 
spect to previous fits, it nevertheless allows a simultaneous estimate 
of R and Ec„,. We found that the reflection parameters of 5 objects 
and the cut-ofi" energy of 4 are well constrained; in the cases of B3 
0309+41 1, 3C 1 1 1 and IGR J21247+5058, these represent the best 
fit results. Due to the above considerations, we employ in the fol- 
lowing discussion the values reported in table 6, with the eventual 
use of the upper liinits on R and £„„ obtained for IGR J13 109-5552 
and S5 2116+81. It is important to note that the values of the pa- 
rameters listed in table 6 are similar to those obtained in previous 
fits and also that in most sources the cross-calibration constants are 
now compatible with 1, an indication that the model is appropriate 
and that variability is not s trong or common at high energies (see 
also lBeckmaim et al]|2007h . 



5 COMPARISON WITH PREVIOUS MEASUREMENTS 

Five out of eight sources in our sample have been previously stud- 
ied over a broad energy range similar to ours, so that a direct com- 
parison is possible. Overall, we find that our parameter values are 
in good agreement, or compatible within errors, with those found 
by these previous works. In the case of 3C 111, all observations 



point to a small reflection and a compatibly low iron line EW 
jGrandi et al. 20Q6; Lewis et al.. 2005); our study yields a slightly 
higher value of R which is however still compatible, within er- 
rors, with previous results. We were able to put a stringent con- 
straint on the high energy cut-olf which we locate at a slightly lower 
value than in previous works, although well within the uncertainties 
found for example by BeppoSAX CGrandi et al..2006 ). 

Comparison o f ou r results for 3C 390 .3 with the works by 
iGrandi et d](ll999') and G liozzi et d]( l2003l) indicates good agree- 
ment on the cut-off energy, although we, as Gliozzi et al., could 
not place a constraint on this parameter, which is instead found by 
Grandi and collaborators to be around 250 keV, i.e. still compat- 
ible with our lower limit of 300 keV. Our value of the reflection 
parameter is also fully compatible with the result obtained using 
BeppoSAX data ( Grandi et al.,. 1999) . but higher (although compat- 
ible wit hin relative uncerta inties) than that obtained with RXTE 
spectra ( iGUozzi et ani2003h . 

We also find very good agreement with previous works for 
4C 74.26 (e.g. iBallantvne & Fabianlbool ICrandi et alj|2006l and 



iLarsson et alj|2008[ ). even though we find a slightly lower cut-off 
energy value. Not much can be said about S5 2116+81 because of 
the poor statistical quality of our data; however, comparison with 
B eppoSAX data shows that our values do not contradict those found 
bv lorandi et"aLH2006l) . 

Finally our results are in full agreement with a previous anal- 
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channel energy (keV) 



channel energy (keV) 



Figure 4. Model to data ratios for 3C 390.3 (left panel) and 4C 74.26 (right panel). The model employed is a simple power law absorbed by both Galactic and 
intrinsic column density (only in the case of 4C 74.26) plus a narrow (3C 390.3) or broad (4C 74.26) Gaussian line component (see Table 3). 



iGRJ21 247+5058 





channel energy (keV) 



channel energy (keV) 



Figure 5. Model to data ratios for S5 21 16+81 (left panel) and IGR J21247+5058 (right panel). The model employed is a simple power law absorbed by both 
Galactic and intrinsic column density plus, in the case of IGR J21247+5058, a narrow Gaussian line component (see Table 3). 



ysis of X-x&y I INTEGRAL data of IGR J21247+5058 jMolina et alj 
l2007i) . For QSO B0241+62, B3 0309+411 and IGR J13109-5552 
no previous work is available and thus the analysis presented here 
is the first broad-band study for those 3 AGN. 



6 DISCUSSION AND CONCLUSIONS 

In the following we discuss the distribution of the fit parameters in 
our sample of AGN (see table 6) in comparison to a similar study 
made on a set of 9 radio quiet Seyfert 1 (Panessa et al. 200^ also 
detected by INTEGRAL. The power law slopes are found to be simi- 
lar in both samples, i.e. around 1.5-2.0; only one of the radio bright 
AGN (the one requiring complex absorption) shows a flat power 
law continuum, while low values of T are more numerous in the 
sample of Panessa et al. ( 2008). The intrinsic absorption measured 
in our sources is generally small or absent (<2xl0^' atoms cm"^), 
except for 3C 111 and IGR J21247+5058, and also similar to the 
values found for radio quiet type 1 AGN; in 3C 111 the extra ab- 
sorption is probably due to intervening material between us and the 
source, while in the case of IGR J2 1247+5058 complex absorption 
is strongly required by the data. It is interesting to note that similar 
complexity, i.e. one or more layers of cold material partially cover- 
in g the central nucleus, is also observed in 2 of the 9 objects studied 
bv lPanessa et al. Il l2008l) . We also point out that, up to very recently, 



the only other BLRG known to have a similar spectra l complexity 
regarding the absorption was 3C 445 dSambruna et al., 2007). for 
which three layers of cold material were required to model the spec- 
trum. Our analysis suggests that objects with complex absorption 
are present within the population of br oad line AGN , indep endent 
of them being radio quiet or radio loud. iGrandi et al.l ( |2006|) found 
evidence in their data for a correlation between absorption and ra- 
dio core dominance, with more absorbed objects showing lower CD 
values. We do not find such a trend; furthermore we do not have 
any indication of a correlation between photon index or absorp- 
tion and Khx using both our sources and the radio quiet Seyferts of 
IPanessa et al] ( |20()&) . 

From the spectral analysis presented here, it is also evident 
that the high energy cut-off (although not required by all the sources 
in the sample) spans a wide range of values from ~40 keV up to 
more th an 300 keV, as alrea dy f ound for radio quiet AGN (see for 
instance iMolina et al. I I2OO6I and IPanessa et a i] |2008l) . We also ex- 
plored the possibility of a correlation between the high energy cut- 
off and the photon index but found nothing; a trend of increasing 
cut-off ener gy with higher F ha s been reported in the literature (see 
for instance fPerola et al]|2002h . but due to the strong dependency 
between these two parameters in the fitting procedure, it is difficult 
to discriminate between any true correlation and induced effects. 

As far as the reflection fraction is concerned, we find in all but 
two cases (namely B3 0309+411 and 4C 74.26) very low values 
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Figure 6. Model to data ratio for 3C 390.3 (left panel) and IGR J21247+5058 (right panel). The model used for 3C 390.3 is a simple power law absorbed by 
Galactic column density plus a blackbody component to model the soft excess and a narrow Gaussian line component (Table 3A). For IGR J21247+5058 the 
model used is instead a simple power law absorbed by Galactic column density and two layers of cold absorbing material partially covering the central source 
plus a narrow Gaussian line component (Table 3B). 



for this component, typically below 1 ; as expected in these sources, 
the EW of the iron line is <100 eV. Indeed if the iron line emis- 
sion is associated with the optically thick material of the disk, one 
would expect the line EW to correlate with the reflection fraction 
as observed. However, in B3 0309+4 11 and 4C 74.26 the observed 
iron line EW is too small for the reflection measured: a possible 
explanation resides in the scatter expected in the EW values due 
to a variation in the iron abundance, or to a possible anisotropy 
of the source s eed photons which might affect the observed spec- 
trum l lPetrucci et al. 2001; Merloni et al. 2006). For a given value 
of R, the EW could also differ according to the value of the 
power law photon index: up to r=2 the EW decreases as the spec- 
trum steepens, while above r=2 the trend is re versed jPerola et al.l 
l2002l :l Mattson. Weaver and & Revnolds|[2007h . B3 0309+411 and 
4C 74.26 both have steep spectra and so a lower EW than expected 
on the basis of the measured R value is not a surprise. 

Overall we can conclude that, as already observed by anum- 
ber of authors (i.e. Grandi et a l. 2006 and Sambruna et al] |l999h . 
the reprocessing features of ttthe radio loud AGN analysed here 
tend to be, on average, quite weak. For a better comparison with 
their radio quiet count erparts, we have ag ain used our data in com- 
bination with those of lPanessa et al.l ('2008). Figure[TT]shows a plot 
of R vs. EW for the entire sample. It is evident that radio loud AGN 
are more confined to a region of the plot characterized by low val- 
ues of EW and to a lesser extent of R, whereas radio quiet objects 
tend to be more spread over both axes. This evidence does not seem 
to be related to the radio core dominance, suggesting that any di- 
lution of the reprocessing features, at least for the sources reported 
here, is not caused by the presence of a jet. Since the observed dif- 
ference is not striking, an accretion flow origin for the X/gamma- 
ray emission is a likely explanation for the production of the re- 
processing features in our radio loud AGN; this agrees well with 
some of them being classified as FR II sources, i.e. the BL RG most 
closely resembling radio quiet AGN tevans et alj|2006h . Within 
this scenario, it is still possible that a different geometry and/or 
accretion flow efficiency involving the disk provides the condition 
for weaker reflection components in BLRG. Alternatively weaker 
reflection features might be the result of reprocessing in an ionise d 
accretion disk, as suggested bv lBallantvne, Ross & Fabi^ ( l2002h ; 
this would alleviate the need for a change in the accretion disk ge- 



ometry and provide a more similar enviroment for both radio loud 
and radio quiet AGN. 

Our sample is of course very small, limited to sources which 
have a small dynamic range of parameters (redshift, loudness, core 
dominance, etc.), and contains two objects (i.e. QSO B0241+62 
and IGR J13109-5552) which deserve further radio studies and are 
not immediately confirmed as BLRG. We hope in the near future 
to improve the statistics, adding newly discovered AGN selected in 
the hard X-ray band in order to verify our findings as well as our 
novel approach of selecting radio loud type 1 sources. 
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Tables 

Spectral Fit Results: wag*wa*Cp<H-zga) 



Name Nh r E„„ a EW Cf C* (doC 

(10^2 cm-^) (keV) (eV) (eV) 



QSO B0241+62 
B3 0309+411 


n 93+0.08 
"■^^-0.08 


'■"--0.07 
1 oT+0.03 

1 m*"'^ 

-Q.Q2 


6.40+°*' 
6 34+8-8? 


lOf 
lOf 


85+« 
100+* 
<30 


1 02+0-22 
-0.18 


1 32+0-32 

4 65+^-^'i' 
+-°^-2.no 


227.1 (234) 
574.6 (603) 


3C111 


Q 43+0.02 
-Q.Q2 


fi 40+0.05 
"•^-0.05 


lOf 


44+0-05 
-0.04 


54+0.11 
"•^ -0.10 


1245.7 (1498) 


IGRJ13109-5552 


-0.27 


170+0.24 
1 74+8S 










2 39+1 84 


41.2(55) 


3C 390.3 


fi44+0.02 


lOf 


71+8 


90+"-'" 


94T8:?8 


2200.4 (2234) 


4C 74.26 


-0.01 


1 7-i+8:8i 


fi 44+8:81 

-0.06 


183!g 


103lg 


54+"-"9 


0.87+|4 


1992.0 (2060) 


S5 2116+81 


<0.19 


1 07+0.12 
'•"-0.12 








1 co+0.60 
1-^^-0.44 


2.76+1-34 


89.1 (128) 


IGRJ21247+5058 


Qgg+0.02 
"•"^-0.02 


1 33+0.02 
'•^^-o.oi 


6 39+005 
-0.08 


lOf 


<30 


32+"-"2 
"■-^ -0.02 


36+^-^^ 
"•-^ -0.02 


2692.3 (2559) 



A: cross-calibration constant between XMM and BAT; B: cross-calibration constant between XMM and IWEGRALflSGRl. 



Table 3A 



3C 390J Spectral Fits Results: waj*Cbb+po+zga) 



(1022^-2) 


r 


kT 
(keV) 


Et. 

Ime 
(keV) 


EW 
(eV) 


^1 


^2 


(dof) 




1 89+"-"2 
'•°'-0.02 


2 41+0.22 


fi 42+0.04 
-0.03 


49+' 


1 37+0.16 
■^'-0.15 


1 5S+0.24 


1970.4 (2232) 



line width fixed to 10 eV. 

A: cross-calibration constant between XMM and BAT; B: cross-calibration constant between XMM and INTEGRAL/ISGRI. 



T^ble 3B 



IGR J21247-H5058 Spectral Fit Results: wag*pcfabs*pcfabs*Cpo+zga) 



r 


n' 
"h 


cfi 


"h 
(10^2 cm-2) 


Cf2 


(keV) 


EW 

(eV) 


Cf 


cf 




1 73+0.05 
'■'^-0.04 


lo.siiJ:** 


n 38+0-03 


1-<15 


n c/;+0-04 
0.86_|j„3 


6 39+"-"' 


<30 


0.6212;« 


7n+0.08 
" -0.06 


2356.2 (2556) 



^ : line width fixed to 10 eV. 

A: cross-calibration constant between XMM and BAT; B: cross-calibration constant between XMM and INTEGRAL/ISGRI. 



Spectral Fit Results: wa^*wa* (cutoff pl+zga) 



Name 


Nh 
(10^2 cm-2) 


r 


EcM 

(keV) 


(keV) 


0- 

(eV) 


EW 
(eV) 


^1 




/ (dof) 


Prob.t 


QSO B0241+61 
B3 0309+411 

3C 1 11 

IGR 113 109-5552 


n 1 s+0.09 

047 +0-02 
"■+--0.02 
<0.46 


1 ^s+0.09 

'■'''-o.o? 
1 on+"-"5 

■ +8:8^ 
'■"■'-o.os 

1 55+"^3' 


>78 

>43 
110+118 
" -40 

>58 


6 40+"-07 
"■+"-0.05 
6 34+"-"7 
"■^+-0.09 

6 40+"-"= 
"■+"-0.05 


lOf 
lOf 

lOf 


qq+46 

"-47 
<30 


1 14+0.23 
-0.21 

56+0.09 
"■■^"-0.08 


1 54+0-39 
•- ^ -0.34 

6 75+3-44 

74+0.I9 
"■'+-0.17 

2 911+1-84 


222.9 (233) 
573.9 (602) 
1234.2 (1497) 

39.8 (54) 


96% 
91% 

99.9% 
83% 


4C 74.26 


0,1 !:;;■;» 


I:™+8f 


>78 


6,44:;;-;;;; 


1 84:;;* 






l.iw;"--!:* 


1989,9 (2(1.S9) 


86% 


,S.-^ 21 K.-SI 


. II.IS 


'■"^:;;.ii 


-Si 












S'),2 i 1271 





A: cross-calibration constant between XMM and BAT; B: cross-calibration constant between XMM and INTEGRAL/ISGRI. 
t : fit improvement with respect to table 3. 

Table 4A 



3C390.3 Spectral Fits Results: wa^*Cbb+cutof^l+zga) 



Nh 

(lO^^cm-^) 


r 


kT 
(keV) 


E„, 
(keV) 


(kiv) 


EW 
(eV) 


^1 


C* 
^2 


A^(doO 


PrBb.t 




1 89+"-"^ 

'■°'-o.oi 


2 41+0.21 
"'■^ -0.16 


NC 


6 43+"-"3 
"■+^-0.04 


48+' 
+°-I4 


1 37+0.16 
'■^'-0.15 


1 ^0+0.23 


1970.2 (2231) 


37% 



^ : line width fixed to 10 eV. 

A: cross-calibration constant between XMM and BAT; B: cross-calibration constant between XMM and INTEGRAL/ISGRI. 
^ : fit improvement with respect to table 3 A. 



Table 4B 



IGR J21247+505S Spectral Fit Results: wag*pcfabs*pcfabs*Ccutof:^l+zga) 



r 


(lO^^cm-^) 


<fi 


(10^2 cm-2) 


d2 


Ec 
(keV) 


4 

line 
(keV) 


EW 
(eV) 




^2 


(<lol) 


Pn)b.t 


1 49+0.05 
■™-0.07 


"■<:!? 




93+1.70 
*-^'-2.05 


Q 27+0-04 
"■'-0.05 


sfi+100 
*"-63 


.39+0.07 
"■^'-0.08 


<30 


Q„+0.05 
"■"-0.06 


n sfi+0.07 
"■'^-0.07 


2276.8 (2555) 


>99.9% 



* line width fixed to 10 eV. 

A: cross-calibration constant between XMM and BAT; B: cross-calibration constant between XMM and INTEGRAL/ISGRI. 
^ : fit improvement with respect to table 3B. 



10 M. Molina et al. 



Table S 

Spectral Fit Results: waj*wa*Cpexrav+zga), 0sR«2, E, =10000 



(10^^ 



(keV) 



(eV) 



EW 
(eV) 



(dof) 



QSOB0241+61 
B3 0309+4 11 
3C 111 

IGRJ13I09-5552 
AC 74.26 
S5 2116+81 



0.48!S;» 
0.26!»;| 

<0.22 



1 g3+0.07 

1 gi+8:J5 

'■^^-0.04 


>0.37 


6 40+"°* 


lOf 


67!*? 


59+"-28 
"■^'-0.11 


>l.20 


^ ^3+0.09 
''■33-p.lo 


lOf 


70+*' 




1 70+0.05 


1 06+'"53 


6 40+"0' 


lOf 


<30 


30+" "'' 
"■^"-0.04 


175+0.23 


>0.1 










I 04+8:8? 


1 30+0 67 
-0.51 


6 45+004 


<186 


51+37 
■"-16 


033+0.07 


2 03+" ™ 


>0.1 








1 10+^-'^ 
'■'"-0.43 



074+0.36 

2 / j+'-i^ 
^■+^-0.97 


224.1 (233) 
560.4 (602) 


92% 
99.99% 


0.35+"" 


1238.5 (1497) 


99.7% 


1 30+--'^'^ 
'■^"-0.74 


40.5 (54) 


66% 


54+0.16 
"■^*-0.17 


1977.0 (2059) 


99.99% 


1 57+151 
'•^'-0.87 


88.5 (127) 


64% 



A: cross-calibration constant between XMM and BAT; B: cross-calibration constant between XMM and INTEGRAL/ISGRl. 
^: fit improvement with respect to table 3. 

Table 5A 



3C 390.3 Spectral Fits Results: wa^, *wa* Cbb+pexrav+zga) , 0^R<2, £^=10000 



(lO^^cm-^) 


r 


kT 
(keV) 


R 


4„e 

(keV) 


EW 
(eV) 


Cf 


^2 


(dof) 


Prob.t 




1 89+0.03 
'■° -0.02 


2.38!g;g 


70+"--"' 
"■' -0.59 


6 43+"03 
°-*^-0.04 


44+10 
-13 


0.88!S;« 


53+0.55 
"■^^-0.22 


1965.4 (2231) 


98% 



- : line width fixed to 10 eV. 

A: cross-calibration constant between XMM and BAT; B: cross-calibration constant between XMM and INTEGRAL/ISGRI. 
' : fit improvement with respect to table 3A. 



Table 5B 

IGR J21247+S058 Spectral Fit Results: wa^*pcfabs*pcfabs*Cpexrav+zga),0<R<2, =10000 



r 




cf, 


(10-- cm"-) 


cf2 


R 


(keV) 


EW 
(eVl 


Cf 




(dof) 


Prob. ' 


, 77+0.05 

'■"-0.05 


8 14+'-" 


36+""3 
"■^"-0.04 


55+0.16 
-0.18 


n 8R+0-05 


93+0-59 
"■^--0.45 


6 -is+0.08 
'■^ -0.08 


<30 




,0+0. 10 
"-"-0.08 


2341.7(2555) 


99.99% 



^ : line width fixed to 10 eV. 

A: cross-calibration constant between XMM and BAT; B: cross-calibration constant between XMM and INTEGRAL/lSGRl. 
' : fit improvement with respect to table 3. 



Table 6 



Spectral Fit Results: wa^^^ *wa* (pexrav+zga) 



Name 


(10-2 ^^-2 


r 


R 


Ec 
(keV) 


EW* 
(eV) 


Cf 


Cf 


X' (dof) 


QSOB0241+61 
B3 0309+4 11 
3C 111 
3C 390.3 
4C 74.26 

IGR J21247+5058 


0,1+0.15 
" -'-0.10 

n d/;+0.03 
0.46_oo3 

14+0.02 
"■"-0.03 
complex* 


1 64+0.22 

1 gg+Ut 
'■ -0.08 
1 73+0.06 

1 89+'^-'^^ 
'■"'^-0.02 
1 70+O.O6 
'■' -0.07 
1 48+006 
'■'°-0.06 


56+"-'*2 
"■•^ -0.41 
3.48+2-'' 

85+i.^^ 

n /;n+8.^§ 
0.60„^ 

1 97+0.69 
'■^^-0.70 

<0.2I 


>86 
35!?' 

>300 
lOO+g" 
79+23 
' -15 


59+42 

<30 
41+12 
*'-ll 
88+23 

<30 


0.85!g;« 

0.40+0" 

95+''-'''^ 
"■'^-0.25 
95+0.17 
"■'^-0.25 

65+""5 
''■°^-0.08 


1 12+0.65 
A CQ+^.^^ 

6.89_35n 

0.52+"-l' 

1 Gl+'i-J^ 
'■i"_0.3n 

1 01+1-39 
0.86+8| 


222.2 (233) 
554.2 (601) 
1227.9 (1497) 
1966.3 (2231) 
1976.0 (2060) 
2276.7 (2555) 



' : line parameters fixed at value obtained in table 3. 
A: cross-calibration constant between XMM and BAT; B: cross-calibration constant between XMM and INTEGRAL/ISGRl 

*N1 =786+2.02 f Q 77+0.04. ist2 =q 77+0.18 cf-)=0 89+0-^0 
■ '■ -1.66' '^^l ^ - -0.05'^^// ^-"-0.13' '^^2-u.ay_Qo^. 



QSO B0241 +62 





Figure 7. Table 6 model for QSO B0241-F62 (left panel) and B3 0309-F411 (right panel). The model is a cut-off power law absorbed both by Galactic and 
intrinsic column density reflected by neutral material plus a narrow Gaussian component describing the iron line. 
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Figure 8. Table 6 model for 3C ill (left panel) and 3C390.3 (right panel). The model is a cut-off power law absorbed both by Galactic and intrinsic column 
density reflected by neutral material plus a narrow Gaussian component describing the iron line. 



4C 74.26 IGRJ21 247+5058 




10 100 ^ 10 100 

Energy (keV) Energy (keV) 



Figure 9. Table 6 model for 4C 74.26 (right panel) and IGR J21247+5058 (left panel). The model is a cut-off power law absorbed both by Galactic and 
intrinsic column density reflected by neutral material plus a narrow Gaussian component describing the iron line. 
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Figure 10. Table 3 model for S5 21 16+81 (left panel) and IGR J13109-5552 (right panel). The model is a simple power law, absorbed both by Galactic and 
intrinsic column densities. 
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Figure 11. Reflection fraction vs. EW for the sample analysed here and the sources presented in lPanessa et al ] 12008^. Blue squares are radio quiet sources, 
while radio loud sources are represented by magenta circles. Arrows represent upper or lower Umits on the parameter values. 
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